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In this work, we extract the QCD Equation of State (EoS) using exper- 
imental results of the 4> meson produced in nuclear collisions at AGS, SPS 
and RHIC energies. The data are confronted to simple thermodynamic ex- 
pectations and lattice results. The experimental data indicate a first order 
phase transition, with a mixed phase stretching energy density between ~ 
1 and 3.2 GeV//m 3 . 

PACS numbers: 25.75.Nq, 12.38.Mh 

1. Introduction 

High energy heavy-ion collisions provide an unique opportunity to test 
the Quantum Chromodynamics (QCD) prediction of a phase transition from 
hadronic matter to a deconfined thermalized state of quarks and gluons 
called Quark Gluon Plasma (QGP). The critical temperature for this tran- 
sition is T c ~ 170 MeV pQ. The order of such a phase transition is still a 
matter of debate. Current lattice QCD calculations indicate that the order 
of the transition depends on the quark masses, as well as on the baryochem- 
ical potential (jib)- The magnitude of the baryochemical potential at the 
central rapidity region depends on the collision energy, it reduces with the 
increase in beam energy. Therefore, by changing the beam energy it is 
possible to study the change in the nature of QCD phase transition. 

In a first order phase transition scenario, the pressure increases with in- 
creasing temperature, until the transition temperature T c is reached, then 
it remains constant during the mixed phase, and continues to increase after 
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the end of the mixed phase. Related to this picture, L. Van Hove [2] sug- 
gested to identify the deconfinement transition in high energy pp collisions, 
looking at the variation of average transverse momentum ((pt)) of hadrons 
as a function of the hadron multiplicity at midrapidity (dN/dy) and search- 
ing for an increase of (pr) followed by a plateau-like behavior and again a 
subsequent increase. The (pt) is expected to reflect the thermal freeze-out 
temperature (Tth) and a flow component which can be related to the ini- 
tial pressure, similarly the hadronic multiplicity reflects the entropy density 
of the system. While the purely thermal component of (pr) can not be 
related to the initial temperature which remains unmeasurable above T c , 
however the flow component in the inverse slope can reflect the plateau of 
the pressure during mixed phase. 

It has been observed that in central heavy ion collisions, the (%) (= 

\jp\ + m 2 ) of pions, kaons and protons as a function of dN^/dy shows a 
Van-Hove-like behaviour as explained above for a wide range of collision 
energies OH]. Hydrodynamic calculations assuming a first order transition 
could reproduce these data [21 [5]. In this work, we study for the first time 
the variation of the inverse slope, T g g extracted from the px distribution of 
the (ft- meson as a function of the initial energy density, e Bj evaluated within 
the framework of Bjorken's hydrodynamical model [B] and ^/sJm for ener- 
gies spanning from AGS, SPS to RHIC [7J. This is an important new feature 
of such studies, because it connects the inverse slope with a parameter char- 
acterizing the initial state of the collision build up after ~ 1 fm/c and which 
reflects at the same time the collision energy, the stopping and the impact 
parameter of the collision. For example, at a given ^/stvtv, different energy 
densities could be achieved by changing the colliding nuclei species or the im- 
pact parameter. €Bj can be directly compared to the critical energy density 
e c obtained in lattice QCD calculations- e c ~ 1 GeV/fm 3 . In the follow- 
ing the (ft data will be analysed and confronted to simple thermodynamic 
expectations which relate to the Van Hove signature and to lattice QCD 
predictions. Due to its ss valence quark content, the (ft meson is of special 
interest [8] to study strangeness enhancement, which is a potential probe of 
QGP formation. Because of its small hadronic rescattering cross section of 
cr(0N) = lO mb, it decouples earlier than other hadrons. Furthermore, due 
to its life time of ~ 45 fm/c, its main decay product (K + K~) suffer less 
rescattering. Experimental results from Au+Au collisions at RHIC energies 
indicate that (ft has a higher thermal freeze out temperature as compared to 
pions, kaons and protons. In particular, their thermal freeze out tempera- 
ture is within errors compatible with the chemical freeze out temperature of 
hadrons and the critical temperature. Therefore (ft and its flow phenomena 
are particularly interesting probes for studying the EoS and the nature of 
the phase transition. 
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V^(GeV) T tfe (MeV) v r 



4.87 150 0.5 

17.3 170 0.5 

62.4 170 0.6 
130 170 0.65 
200 170 0.65 



Table 1. Freeze-out temperature, Tth and the radial flow velovity, v r for 0- meson 
for various centre of mass energies, ^snn extracted from the experimental data 
within the framework of blast wave method. 
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Fig. 1. Transverse mass distribution of as a function of its kinetic energy at 
mid-rapidity for different colliding energies measured by the STAR collaboration 
at RHIC (left) and for SPS and AGS energies (right). The theoretical curves are 
obtained within the framework of blast wave method. 



2. as a probe of the order of the phase transition 

The pt spectra of <j) measured in heavy ion collisions at various collision 
energies [9j fl"0| [TT] have been reproduced within the ambit of the blastwave 
method [16], results for AGS, SPS and RHIC are shown in figure [TJ The 
values of the parameters i.e. the radial flow velocity, v r and the freeze- 
out temperature, of the blastwave method are displayed in table I (see 
also |18]). It has been found that Tth ~ 170 MeV and v r ~ 0.6 and are 
independent of the center of mass energies beyond a thershold. T t h being 
close to T c indicates that <f> meson freezes out near the phase boundary and 
could be used to extract the properties of QCD matter near the transition 
point. Moreover the values of v r indicate that QGP has undergone substan- 
tial radial flow. The values of the 'true' temperature and the flow velocity 
of the system at the <f> freeze-out surface for different colliding energies can 
be used to get the effective inverse slope parameter T e ff = T t h + \m (v r ) 2 . 



4 



excitedQcdAPPBFI printed on August 21, 2009 



450 

— 400 
a 

— 350 
300 
250 
200 
150 




^] (GeV) 1 



Fig. 2. Tf ff (left ) and (right) as a function of </sWn from AGS-SPS to RHIC. 
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Fig. 3. Tf ff and as a function of egj.r from AGS-SPS to RHIC. 



In figure [2] the effective temperature of <fi, T^g- (left) and the thermal 

component of the inverse slope (right) are shown as a function of the 
collision energy ^/snn- We use the central collision data in almost the same 
Pt range, at mid-rapidity [91 [TO], [Til 021 13 O [H] . It is observed that from 
AGS to SPS energies these observables remain almost unchanged, showing a 
plateau-like structure. Going from SPS to RHIC energies, the inverse slope 
T'eff exhibits a sudden jump while an increase is still observed in the 
component. This may be due to an imperfect transverse flow component 
subtraction at RHIC or other effects. It is observed that the is reaching 
at RHIC values compatible within errors with T c . The observed plateau 
of T^q is a signature of a coexisting phase of quarks, gluons and hadrons 
for a first order phase transition, during which the initial pressure remains 
constant. The subsequent increase of with ^Jsnn at top RHIC energies 
indicates the end of the mixed phase and the entering into a pure QGP 
phase. 

Now we study the dependence of the inverse slope and its thermal com- 
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ponent on the initial energy density, e Bj , estimated as e Bj = 




tkK 1 



1 



( m T) 7^ where R = RqA 1 / 3 and A ~ N part /2. From the exper- 
imental measurements of dE^jdy and dN/dy with (tut), the observable 
e_Bj could be estimated for all centralities and center of mass energies at 
mid-rapidity. 

In figure [3] (left) and (right) are shown as a function of esj-T. 
Note that the formation time, r is model-dependent and in general, in sub- 
sequent discussions we also assume r ~ 1 fm/c. The above result reflects 
the properties of equation of state. We observe a plateau in T^q stretch- 
ing between esj ~ 1 and 3.2 GeV//m 3 , and increasing suddently above 
3.2 GeV//m 3 . This behaviour as already discussed, suggests a first order 
transition, however from figure Owe can now infer that the mixed phase is 
stretching between 1 and 3.2 GeV/fm 3 . The use of the esj scale, allows 
us to establish here for the first time the esj range of the mixed phase. 
More data on the <j) at esj below 1 GeV//m are needed to establish the 
increase of T^jj up to 1 GeV//m 3 , as seen in other hadrons. The increase 
of the inverse slope at RHIC energies again indicate a pure phase of QGP, 
as is expected from a first-order phase transition. The thermal component 
"T'tn shows also a plateau while a smaller increase is still observed at 3.2 

GeV/ fm 3 . The step-like behavior in the excitation function of < mrp > 
has also been explained by taking a first-order phase transition with a large 
latent heat or if the EoS is effectively softened due to non-equilibrium effects 
in the hadronic transport calculations |17j . 

The flow component of the inverse slope of the 4> reflecting the initial 
pressure, is of interest to look directly observables linked to this initial 
pressure like the transverse flow velocity, v r for all hadrons and the elliptic 
flow t>2 a s a function of collision energy. The transverse flow velocity v r 
shows exactly the same characteristic behaviour as the T^g. as a function 
of collision energy, namely increase, a plateau and subsequent increase at 
RHIC |18] . A similar pattern is suggested for the elliptic flow V2 as a function 
of collision energy [TS]. The observation of a mixed phase goes inline with 
the fact that excitation function of various observables show anomalous 
behavior or saturation effects starting lower SPS energies \20\ I21j . This 
could be related to the onset of deconfinement corresponding to this energy 
regime. In addition, the fact that microscopic transport models, based on 
hadronic degrees of freedom failed to reproduce the observed behavior of 
kaon inverse slope \22\ [23] . also indirectly confirms the observation of a 
deconfinement transition. 
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Fig. 4. Energy density as a function of Tth (left) and T e ff (right) for the <fi meson. 



3. Comparison of data with lattice predictions 

In the following, we compare experimental data with lattice QCD predic- 
tions on energy density as a function of temperature as shown in figure 8 
of [21]. In figure [U the lattice prediction for the energy density as a func- 
tion of the temperature is compared to eBj as a function of Tf h (left) and 
Tfjrjr (right) extracted from the slopes of <fi spectra. First we discuss the left 
plot. The lattice prediction is for zero baryochemical potential therefore 
corresponds to the cross over region. The figure in the left is at non-zero 
[iB and shows the variation of the initial energy density with temper- 
ature, which is measured at a later time than the energy density, namely 
at the thermal freeze out of the 4>. That is, the two variables in the data 
are a measure of the system at different times, while the lattice estimate is 
independent of time. This temperature is expected to be always below T c 
and does not reflect the initial T, which may exceed T c depending on the 
collision energy. This plot can be directly compared to figure 11 of [25], 
where the energy density from data has been studied as a function of T, 
while both x and y axis were estimated at the same time, namely at the 
hadronic chemical freeze out time and at = 0. It is seen that the energy 
density increases approaching the T c from below, from which both the T c 
and critical exponents can be extracted [25] . 

If the transition occurs at T c ~ 170 MeV as expected, the temperature 
of <j> at the thermal freeze out should saturate below and near T c , for all 
values of the initial density (up to infinity) , as seen in figure U] (left panel) . 
We observe a saturation in the value of Tf h near T c , because Tj h < T c as 
mentioned earlier. To measure temperatures above T c photons and dileptons 
are useful probes. 
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We now discuss figured] (right panel), where the initial Bjorken energy 
density is shown as a function of the effective temperature of the <f> for 
collisions from AGS, SPS to RHIC. Now both the x and y axis are reflecting 
parameters at initial times therefore this plot is more appropriate to be 
compared to the lattice results. The effective temperature of the (f> here is 
a sum of a thermal freeze out temperature, which is expected to be always 
below T c , and a non-thermal component due to transverse flow, which relates 
to the initial pressure and reflects the initial conditions above T c . Therefore 
the variables in the two plots compared here are not exactly the same but 
they are correlated. Further analysis is needed to compare the data to 
lattice e.g. using exactly the same variables in both estimates. A study can 
be also done involving the elliptic flow V2 and the transverse flow velocity v r 
as a function of the initial energy density. Also a hydrodynamic calculation 
of the discussed variables is of interest. 



4. Summary 

In summary, a first analysis of experimental data on the inverse slope 
parameter of <p mesons as a function of ^Jsnn & energy density from AGS, 
SPS to RHIC, and their comparison to simple thermodynamic expectations, 
suggest a first order phase transition as predicted by lattice QCD at non-zero 
baryochemical potentials. The mixed phase of quarks, gluons and hadrons, 
is found to stretch between esj of ~ 1 and 3.2 GeV//m 3 . The latent heat 
density in a first order phase transition is eg(T c ) — e#(T c ) = 4B, where 
€q(€h) is the energy density of QGP (hadrons) at T c and B is bag constant. 
With e Q {T c ) ~ 3.2 GeV/fm 3 and e H {T c ) ~ 1 GeV/fm 3 , we get a reasonable 
value for B 1 ^ ~ 250 MeV. If we ignore B, then €q/€h ~ 9q/9h ~ 3.2 
also a reasonable number comparable to results from lattice QCD [T], here 
9q{9h) is statistical degeneracy of QGP (hadrons). The plateau and subse- 
quent increase of the inverse slope parameter of cj) above eBj 3.2 GeV/fm 3 
is in agreement with data from pions, kaons and protons, and is as well 
observed in the transverse flow velocity v r reflecting the behaviour of the 
initial pressure. A first attempt to compare data to lattice QCD predictions 
is made. 

The above results, while supporting the order of the transition predicted 
by lattice at large \xb up to esj ~ 3.2 GeV/fm 3 , they do not exclude a 
change of order of the transition at smaller baryochemical pontials. Which 
points to further work, towards mapping out the order of the QCD phase 
transition as a function of energy and baryochemical potential. 
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